Thromboembolic disease markers

ABSTRACT

The invention relates to a method for a more appropriate thromboembolic event risk assessment based on the presence of different genetic variant. The invention also relates to a method for determining the risk of suffering a thromboembolism disease by combining the absence or presence of one or more polymorphic markers in a sample from the subject with conventional risk factors for thromboembolism as well as computer-implemented means for carrying out said method.

RELATED APPLICATIONS

This application is a continuation-in-part of U.S. application Ser. No. 14/886,463, which is a divisional of U.S. application Ser. No. 14/126,624, filed Mar. 4, 2014, which is a national stage filing under 35 U.S.C. § 371 of international application PCT/EP2012/061185, filed Jun. 13, 2012, which was published under PCT Article 21(2) in English, the disclosure of each of which is incorporated by reference herein in its entirety.

FIELD OF THE INVENTION

The present invention relates to the field of thromboembolic diseases or disorders. More specifically, it relates to markers and methods for determining whether a subject, particularly a human subject, is at risk of developing thromboembolic disease or disorder, developing a thromboembolic event, having thromboembolic disease or disorder, or experiencing a complication of a thromboembolic disease.

TECHNICAL BACKGROUND

Thromboembolic disease is the leading cause of morbidity and mortality in the developed world (America Heart Association 2010. Circulation 2010; 121:e46-e215). Arterial thrombosis is the most common underlying cause of acute myocardial infarction, non-hemorrhagic cerebrovascular accidents, and peripheral vascular disease. Pathological manifestations of venous thromboembolism (VTE) largely include deep vein thrombosis (DVT) and pulmonary embolus (PE). While arterial thromboembolic events are the foremost cause of death and disability, venous disease also plays an important role. VTE occurs for the first time in 100 per 100 000 persons each year in the United States (America Heart Association 2010. Circulation 2010; 121:e46-e215). Approximately one third of patients with symptomatic VTE manifest PE, whereas two thirds manifest DVT alone (America Heart Association 2010. Circulation 2010; 121:e46-e215). PE is the most common cause of preventable hospital death accounting for 60.000 deaths in the United States annually (Anderson F A. Arch Intern Med 1991: 151:933-8, Spencer F A. Arch Inter Med 168:425-430).

Medical textbooks and epidemiological studies characteristically consider arterial and venous thromboembolic disease as distinct entities, each with their own pathophysiological basis, unique risk factors, and distinct therapeutic regimens (Bauer K A. Hematology Am Soc Hematol Educ Program 2002; 353-368). Arterial clots typically occur in an injured vessel and the most common cause of vascular damage in the arterial system is atherosclerotic vascular disease (AVD) (Lane D E 2000. Thromb Haemost 2000; 76:651-62). The risk factors for arterial thrombosis are therefore considered the same as those for AVD. Arterial clots occur in a high flow, high shear environment and these clots, also called white clots, are rich in platelets. Prevention and treatment of arterial thrombosis is often aimed at platelet inhibition. While vascular injury can promote the formation of venous clots, stasis and changes in blood composition (thrombophilia) are the most important risk factors for venous clot development (Lane D E 2000. Thromb Haemost 2000; 76:651-62). Venous clots occur in a low flow system, they are rich in fibrin that is enmeshed with red blood cells and are referred to as red clots. Inhibition of fibrin formation is the mainstay of prevention and treatment of venous thrombosis. It is often reported that the risk factors for arterial and venous thrombosis largely differ (Bauer K A. Hematology Am Soc Hematol Educ Program 2002; 353-368).

However, recent studies have demonstrated a close association between arterial and venous thrombosis at a variety of levels. Specifically it has been shown that:

-   1) arterial and venous thrombosis share common risk factors (Doggen     C J M Arteroscler Thromb Vasc Biol 2004; 24:1970-5., Goldhaber S Z     In: Bloom A L, Forbes C D, Thomas D P, Tuddenham E G D, eds.     Hemostasis and Thrombosis. New York: Churchill and Livingstone:     1997; 1327-1333.), -   2) individuals who suffer idiopathic venous thromboembolism are at a     markedly increased risk of suffering a significant cardiovascular     event (Becattini C. European Heart Journal 2005; 26:77-83.), -   3) individuals who suffer idiopathic venous thromboembolism have an     increased incidence of atherosclerotic vascular disease     (Becattini C. European Heart Journal 2005; 26:77-83.), and -   4) those who suffer idiopathic venous thromboembolism have a     significantly higher incidence of metabolic syndrome (Ageno W. J     Thromb Haemost 2006; 4:1914-8).

The risk factors reported to be common to both arterial and venous thrombosis and that represent significant hazard for the development of each entity include increasing age and weight, smoking, exposure to estrogens, and the presence of diabetes. It has also been shown that high HDL cholesterol levels are associated with a decreased risk of venous thrombosis while elevated triglyceride and/or total cholesterol levels convey an increased risk. Other risk factors reported to be common to both arterial and venous thrombosis include the presence of antiphospholipid antibodies, dysfibrinogenemia, hyperhomocysteinemia, and elevated levels of fibrinogen, lipoprotein (a) and factor VIII.

One of the strongest pieces of evidence in favour of a link between arterial and venous thrombosis is the Genetic Analysis of Idiopathic Thrombophilia (GAIT) Study (Souto J C. Am J Hum Genet 2000; 67:1452-1459). This family-based study of the genetics of thrombosis in a Spanish population was initiated to determine the heritability of thrombosis. Three hundred and twenty-eight individuals in 21 extended pedigrees were evaluated using a novel computer assisted adaptation of a multivariate threshold model. The authors concluded that more than 60% of the variation in susceptibility to common thrombosis is attributable to genetic factors. What makes this study unusual is that both venous and arterial thromboembolic events were included in the analysis. When venous and arterial thrombosis were jointly analyzed, arterial and venous thromboses were highly genetically correlated. That is, many of the same genes are involved in the pathogenesis of arterial and venous disease.

There are also studies (Doggen C J M, Smith N L, Lamahre R N, et al. Arteros Thromb Vasc Biol 2004; 24:1970-5; Becattini E. European Heart Journal 2005; 26:77-83; Ageno W.] Thromb Haemost 4:1914-8) suggesting that arterial and venous thrombosis represent different manifestations of the same disease and that the underlying process is driven by a common set of genes.

A) Prevention of First Episode of Thromboembolism

Symptomatic thrombosis (arterial or venous) is a multifactorial disease that manifests when a person with an underlying predisposition to thrombosis (thrombophilia also referred to as thrombophilic disorder or hypercoagulable syndromes) is exposed to clinical risk factors.

Assessment of presence of thrombophilia is not solely confined to laboratory testing but begins with a detailed history and physical examination. Detailed inquiry into symptoms and signs of acquired risk factors (coexisting diseases, medication exposure, and clinical circumstances) that are associated with thrombosis are an important part of the initial evaluation as is a complete physical examination. In addition to judicious laboratory testing appropriate for the patient's age and symptoms, objective confirmation of venous thromboembolism is critical.

Laboratory Testing

Currently, there is no single laboratory global assay that will ‘screen’ for the presence of thrombophilia. Thus, laboratory testing can be broadly categorized into (1) general diagnostic testing, (2) specialized coagulation testing, and (3) ancillary testing for disorders known to predispose to thrombotic disorders.

Specialized Coagulation Testing

Special coagulation testing consists of a battery of complex (protein and DNA-based) thrombophilia assays to detect presence of an inherited or acquired thrombophilia. However, multiple preanalytical conditions affect results of the non-DNA-based assays (e.g. anticoagulants, acute thrombosis, liver disease, etc.), so interpretation of results needs to be done within the context of the circumstances surrounding testing. An additional factor affecting the yield of testing is the ethnicity of the patient population being studied. Prevalence of factor V Leiden (FVL) varies from 3% to 7% in Caucasians of European ancestry, but has a very low prevalence in individuals of other ethnic groups: 0% among Native Americans/Australians and Africans, 0.16% among the Chinese, and 0.6% among individuals from Asia Minar (India, Pakistan, Sri Lanka). The ethnicity is important especially when few (one or two) genetic markers are analysed.

Factors Affecting Results of Protein-Based Specialized Coagulation Testing.

Effect of Acute Thrombosis

During the acute thrombotic episode, levels of antithrombin, protein C, and protein S may be transiently reduced; thus, if testing is not repeated, remote from the thrombotic event and from anticoagulant therapy, the patient may be misdiagnosed as having a congenital deficiency.

Effect of Anticoagulants

-   -   Heparin. Heparin therapy can falsely reduce antithrombin levels.         Although most lupus anticoagulant (LAC) reagents [e.g. dilute         russel viper venom time (DRVVT) and Stadot APTT] contain heparin         neutralizers that can neutralize up to 1 U/mL of heparin,         presence of excess heparin may result in a false-positive test         result, which impacts the duration of secondary prophylaxis.         Thus, positive results of LAC testing performed while on heparin         should be reconfirmed when the patient is off heparin.     -   Vitamin K antagonist (VKA) therapy. Protein C and S levels are         lowered by VKA therapy (e.g. warfarin since they are vitamin         K-dependent proteins). In addition. VKA therapy may result in a         false-positive LAC with certain assays (e.g. DRVVT).     -   Direct thrombin inhibitors (DTIs; e.g. argatroban, lepirudin,         desirudin, bivalirudin, Dabigatran etexilate or Dabigatran).         Because the majority of anticoagulant activity assays rely on         generation of thrombin to achieve an endpoint of dot detection,         presence of DTIs interfere with this endpoint and delay dot         formation. This can lead to a false-positive LAC or falsely         reduced protein C and S levels. Results of chromogenic assays         are likely reliable.     -   Factor Xa inhibitor. Examples of factor Xa inhibitors are         fondaparinux, rivaroxiban, betrixaban, edoxaban, otamixaban,         letaxaban, eribaxaban and apixaban.         Effect of Liver Disease

The majority of anticoagulant and procoagulant proteins are produced in the liver. In advanced liver disease, levels of both the anticoagulant and procoagulant proteins are reduced.

Sample Collection and Processing Issues

Practically speaking, ordering physicians have limited impact on specimen collection and processing; however, knowledge of such effects may lead one to consider repeat testing, if the data are unexpected or do not fit the expected pattern [e.g. reduced activated protein C resistance (APC-R) ratio suggesting presence of APC-R, yet the FVL test is negative].

Effect of Type of Anticoagulant in Specimen Collection Tube

Standard specimen collection tubes contain 0.105-0.109 mol citrate for optimal results. Specimens may inadvertently be collected in ethylenediaminetetraacetic acid (EDTA), which will result in falsely reduced protein levels and a reduced APC-R ratio.

Effect of Specimen Processing

Specimens should be double centrifuged as soon as possible after collection in order to reduce the amount of residual platelets to a minimum. The presence of residual platelets can result in a false-negative test for LAC.

Molecular Risk for Thrombotic Disease

Although an inherited tendency for excessive bleeding is often be ascribed to single or few gene abnormalities, there is ample evidence to suggest that, in contrast, the clinical manifestations of hypercoaguability are usually the result of adverse interactions between multiple genes and the environment. Thus, the use of molecular diagnostics to document markers of thrombotic risk (thrombophilia) will prove to be far more challenging than with the inherited hemorrhagic disorders. To further complicate matters, despite the fact that with appropriate testing, thrombophilic mutations can be identified in patients following a first clinical episode of venous thromboembolism, interpretation of these results remains problematic.

Inherited Resistance to Activated Protein C: Factor V Leiden

Until 1994, the investigation of patients with clinical evidence of hypercoaguability was usually unproductive. However, with the discovery by Dahlback and Hildebrand of an inherited form of resistance to the proteolytic effects of activated protein C, and the subsequent finding of a common missense mutation in the factor V gene by Bertina and colleagues in Leiden, a major advance was made in the laboratory assessment of thrombotic risk.

The Leiden mutation substitutes a glutamine for an arginine at amino acid residue 506 in factor V, the initial cleavage site for activated protein C. The mutation is readily detected by a number of PCR-based approaches. Between 2% and 5% of individuals in Western populations have been documented to be heterozygous for factor V Leiden. In contrast, the mutation is extremely rare in subjects of Asian and African descent.

In some laboratories, initial screening for resistance to activated protein C is performed using the prolongation of an activated partial thromboplastin time-based assay as an indicator; patients testing positive (prolongation in the presence of factor V-deficient plasma) are subsequently evaluated by a PCR.

Increasingly, where access to PCR-based molecular analysis is routine, laboratories will more often choose to proceed directly to the genetic test, as the result is definitive and more than 95% of activated protein C resistance is a result of this single mutation.

Persons heterozygous for the factor V Leiden mutation have an approximately five-fold increased relative risk of venous thrombosis. It is found in 15-20% of patients experiencing their first episode of venous thrombosis. The hypercoagulable phenotype associated with factor V Leiden shows incomplete penetrance, and some individuals may never manifest a clinical thrombotic event. In contrast to the increased relative risk for an initial venous thrombotic event associated with factor V Leiden, this genetic variant is not associated with increased risks for either arterial thrombosis or a recurrence of venous thrombosis. Coinheritance of other inherited thrombotic risk factors or exposure to environmental risk factors can dramatically enhance the thrombotic risk in carriers of factor V Leiden. Many clinicians test for this disorder in patients with a family history of thrombosis who are about to be exposed to an acquired thrombotic risk factor. Individuals homozygous for the mutation have a 70-fold enhanced relative risk of venous thrombosis, indicating that this phenotype is transmitted as a codominant trait.

Prothrombin 20210 3′ Non-Coding Sequence Variant

In 1996, Poort and colleagues described an association between a G to A nucleotide polymorphism at position 20210 in the 3′ untranslated region (UTR) of the prothrombin gene, increased plasma levels of prothrombin, and an enhanced risk for venous thrombosis. This polymorphic nucleotide substitution is at the very end of the 3′ UTR and exerts its effect on prothrombin levels in the heterozygous state. Although the plasma levels of prothrombin in subjects heterozygous for this polymorphism are higher on average than those in individuals with a normal prothrombin genotype, levels are usually still within the normal range. As a consequence, this polymorphism can only be evaluated by genetic testing, which is achieved by a PCR-based assay, most often now involving a form of real-time quantitative assay.

As with the factor V Leiden genotype, the prevalence of the prothrombin 20210 G to A variant in the general population is relatively high at 1-5%. This variant is also rare in persons of Asian and African descent. The heterozygous state is associated with a two- to four-fold increase in the relative risk for venous thrombosis. There is no influence on venous thrombotic recurrence. The relationship of prothrombin G2010A with arterial thrombosis is very modest (OR 1.32; 95% CI 1.03-1.69) (Kim R J. Am Heart J 2003; 146:948-957).

Thermolabile C671T 5,10-Methylene-Tetrahydrofolate Reductase Variant

The third, high-prevalence gene tlc variant that was initially thought to be associated with an increased thrombotic risk is the C to T variant at nucleotide 677 (an alanine to valine substitution) in the 5,10-methylene-tetrahydrofolate reductase (MTHFR) gene. This genotype results in expression of an enzyme with increased thermolability. Homozygosity for the variant is associated with hyperhomocysteinemia, particularly in the presence of folate deficiency. In many populations (southern Europeans and Hispanic Americans), approximately 10% of subjects are homozygous for the C677T variant, a sequence change that can easily be detected by a PCR-based strategy. However, after further extended analysis, in contrast to the factor V Leiden and prothrombin 20210 variants, the role of the MTHFR C677T polymorphism as an independent risk factor for venous thromboembolism appears minor.

B) Diagnosis of Venous Thromboembolism

Objective testing for deep vein thrombosis and pulmonary embolism is essential because clinical assessment alone is unreliable. Failure to diagnose venous thromboembolism is associated with a high mortality, whereas inappropriate anticoagulation can lead to serious complications, including fatal haemorrhage.

Diagnosis of Deep Vein Thrombosis

The clinical features of deep vein thrombosis include localized swelling, erythema, tenderness, and distal edema. However, these features are nonspecific, and approximately 85% of ambulatory patients with suspected deep vein thrombosis will have another cause for their symptoms. The differential diagnosis for deep vein thrombosis includes:

-   -   cellulites;     -   ruptured Baker cyst;     -   muscle tear, muscle cramps, muscle hematoma;     -   external venous compression;     -   superficial thrombophlebitis; and     -   post-thrombotic syndrome.         Venography

Venography is the reference standard test for the diagnosis of deep vein thrombosis. It has advantages over other tests in that it is capable of detecting both proximal vein thrombosis and isolated calf vein thrombosis. However, the disadvantages are that it:

-   -   is invasive, expensive, and requires technical expertise; and     -   exposes patients to the risks associated with contrast media,         including the potential for an allergic reaction or renal         impairment.

For these reasons, noninvasive tests such as venous ultrasonography and D-dimer testing, alone or in combination with clinical assessment, have largely replaced venography.

Compression Venous Ultrasonography

This is the noninvasive method of choice for diagnosing DVT. The common femoral vein, superficial femoral vein, popliteal vein, and proximal deep calf veins are imaged in real time and compressed with the transducer probe. Inability to compress the vein fully is diagnostic of venous thrombosis. Venous ultrasonography is highly accurate for the detection of proximal vein thrombosis with a sensitivity of approximately 97%, specificity of approximately 94%, and negative predictive value of approximately 98% in symptomatic patients. If DVT cannot be excluded by a normal proximal venous ultrasound in combination with other results (e.g. low clinical probability or normal D-dimer), a follow-up ultrasound is performed after 1 week to check for extending calf vein thrombosis (present in approximately 2% of patients). If the second ultrasound is normal, the risk of symptomatic VTE during the next 6 months is less than 2%.

The accuracy of venous ultrasonography is substantially lower if its findings are discordant with the clinical assessment and or if abnormalities are confined to short segments of the deep veins. Ideally, these patients should have a venogram because the result of the venogram will differ from the venous ultrasound in approximately 25% of these cases. If venography is not available, additional testing (e.g. D-dimer, serial venous ultrasonography) may help to clarify the diagnosis and avoid inappropriate anticoagulant therapy.

Venous ultrasonography of the calf veins is more difficult to perform (e.g. sensitivity 70%), and its value is controversial. Some investigators have proposed that a single complete compression ultrasound that includes examination of the calf veins should be used to exclude DVT. Studies using this method have reported an incidence of VTE of 0.5% during 3 months follow-up after a negative examination, establishing that a negative venous ultrasound that includes the calf veins excludes VTE [8]. However, this method has the potential to diagnose calf DVT that would have spontaneously lysed without treatment and to yield false-positive results, thereby exposing patients to the risk of bleeding due to anticoagulant therapy without clear benefit.

D-Dimer Blood Testing

D-dimer is formed when cross-linked fibrin is broken down by plasmin, and levels are usually elevated with DVT and/or PE. Normal levels can help to exclude VTE, but elevated D-dimer levels are non-specific and have low positive predictive value. D-dimer assays differ markedly in their diagnostic properties for VTE. A normal result with a very sensitive D-dimer assay (i.e. sensitivity of approximately 98%) excludes VTE on its own [i.e. it has a high negative predictive value (NPV). However, very sensitive D-dimer tests have low specificity (approximately 40%), which limits their use because of high false positive rates. In order to exclude DVT and/or PE, a normal result with a less sensitive D-dimer assay (i. e. approximately 85%) needs to be combined with either a low clinical probability or another objective test that has a high NPV, but is non-diagnostic on its own (e.g., negative venous ultrasound of the proximal veins. As less sensitive D-dimer assays are more specific (approximately 70%), they yield fewer false-positive results.

Specificity of D-dimer decreases with aging and with comorbid/illness, such as cancer. Consequently, D-dimer testing may have limited value as a diagnostic test for VTE in hospitalized patients (more false positive results) and is unhelpful in the early postoperative period.

Computed Tomographic (CT) Venography and Magnetic Resonance (MR) Venography

CT venography and MR venography have the potential to diagnose DVT in settings where the accuracy oi compression ultrasonography is limited (e.g. isolated pelvic DVT, asymptomatic patients). The sensitivity and specificity of CT venography compared with compression ultrasonography for detecting all DVT has been reported between 89% and 100%, and 94% and 100%, respectively. However, given the cost, exposure to radiation, and limited availability of CT venography, this modality currently plays a limited role in the diagnosis of DVT. A meta-analysis of studies comparing MR venography with conventional venography reported a pooled sensitivity of 92% and specificity of 95% of MR venography for proximal DVT. As with CT venography, cost and availability will inhibit the widespread use of MR for diagnosis of acute DVT.

Diagnosis of Pulmonary Embolism (PE)

The clinical features of PE may include:

-   -   pleuritic chest pain,     -   shortness of breath,     -   syncope,     -   hemoptysis, and     -   palpitations.

As with DVT, these features are non-specific, and objective testing must be performed to confirm or exclude the diagnosis of PE.

Pulmonary Angiography

This is the reference standard test for the diagnosis of PE. However, it has many of the same limitations as venography.

Computed Tomographic Pulmonary Anglography (CTPA)

Spiral CT (also known as helical CT) with peripheral injection of radiographic contrast (CTPA) is the current standard diagnostic test for PE (Stein P D. N Engl J Med 2006; 354:2317-2327, Roy P M. Br Med J 2005; 331:259). In comparison with ventilation-perfusion lung scanning, CTPA is less likely to be “non-diagnostic” (i.e. approximately 10% vs. 60%) and has the potential to identify an alternative etiology for the patient's symptoms. This technique has a sensitivity of 83%, specificity of 96%, NPV of 95%, and positive predictive value of 86% for PE.

Accuracy of CTP A varies according to the size of the largest pulmonary artery involved and according to clinical pretest probability. For example, the positive predictive value of CTPA is 97% for pulmonary emboli in the main or lobar artery, but drops to 68% for segmental arteries, and is lower still for PE in the subsegmental arteries (25%). In patients with a high clinical pretest probability of PE, the positive predictive value of CTPA is 96%, but this value falls to 92% in patients with a technical pretest probability of PE, and to 58% in patients with a low clinical pretest probability of PE.

In management studies that used CTP A to diagnose PE, less than 2% of patients who had anticoagulant therapy withheld based on a negative CTPA went on to have symptomatic VTE during follow-up. Taken together, these observations suggest the following:

-   -   A good-quality, normal CTPA excludes PE if clinical suspicion is         low or moderate.     -   Lobar or larger pulmonary artery intraluminal defects are         generally diagnostic for PE.     -   Segmental pulmonary artery intraluminal defects are generally         diagnostic for PE if clinical suspicion is moderate or high, but         should be considered non-diagnostic if suspicion is low or there         are discordant findings (e.g. negative D-dimer).     -   Subsegmental pulmonary artery intraluminal defects are         nondiagnostic, and patients with such findings require further         testing.

A note of caution: If possible, CTPA should be avoided in younger women (e.g. younger than 40 years) because it delivers a substantial dose of radiation to the chest, which increases the risk of breast cancer.

Ventilation-Perfusion Lung Scanning

In the past, ventilation-perfusion lung scanning was the initial investigation in patients with suspected PE, and it is still useful in patients with contraindications to x-ray contrast dye (e.g. renal failure) and patients at higher risk for developing breast cancer from radiation exposure (e.g. young women). A normal perfusion scan excludes PE, but is only found in a minority of patients (10-40%). Perfusion defects are non-specific; only approximately one-third of patients with perfusion defects have PE. The probability that a perfusion defect is caused by PE increases with size and number and the presence of a normal ventilation scan (“mismatched” defect). A lung scan with mismatched segmental or larger perfusion defects is termed “high-probability.” A single mismatched defect is associated with a prevalence of PE of approximately 80%. Three or more mismatched defects are associated with a prevalence of PE of approximately 90%. Lung scan findings are highly age-dependent, with a relatively high proportion of normal scans and a low proportion of non-diagnostic scans in younger patients. A high frequency of normal lung scans is also seen in pregnant patients who are investigated for PE.

Clinical Assessment:

As with suspected DVT, clinical assessment is useful for categorizing probability of PE.

D-Dimer Testing:

As previously discussed when considering the diagnosis of DVT, a normal D-dimer result, alone or in combination with another negative test, can be used to exclude PE.

Patients with Nondiagnostic Combinations of Noninvasive Tests for PE

Patients with non-diagnostic test results for PE at presentation have a prevalence of PE of approximately 20%; therefore, further investigations to exclude PE are required.

Diagnosis of PE in Pregnancy

Pregnant patients with suspected PE can be managed similarly to non-pregnant patients, with the following modifications:

-   -   Venous ultrasound of the legs should be performed first followed         by ventilation-perfusion lung scanning if there is no DVT.     -   The amount of radioisotope used for the perfusion scan should be         reduced and the duration of scanning extended.     -   If pulmonary angiography is performed, the brachial approach         with abdominal screening is preferred.     -   The use of CTPA in pregnancy is discouraged, primarily because         of radiation exposure to the mother.         C) Risk of Recurrence after a First Episode of Symptomatic         Venous Thromboembolism

Venous thromboembolism is associated with diverse risk factors, some of which are transient, such as recent surgery and pregnancy, and others of which are persistent, such as cancer (Table 1 shows the risk factors for venous thromboembolism). When venous thromboembolism is associated with an acquired risk factor, either transient or persistent, it is called provoked. When there is no apparent clinical risk factor, it is called unprovoked or idiopathic.

TABLE 1 Risk factors for venous thromboembolism Major transient risk factors Hospitalization Plaster cast immobilization Surgery Trauma Minor transient risk factors Oral contraceptives or hormone therapy Pregnancy Presence of major risk factors 1 to 3 months before venous thromboembolism Prolonged travel (≥2 hours) Potential acquired or persistent risk factors Collagen vascular diseases Heart failure Malignancy Medications Myelo proliferative disorders Nephrotic syndrome

It has recently been recognized that the presence or absence of a transient, or reversible, risk factor at the time of venous thromboembolism strongly affects the risk of recurrence after anticoagulant therapy is stopped. Patients with venous thromboembolism provoked by a transient risk factor have a low risk of recurrence compared with patients with either venous thromboembolism provoked by a persistent risk factor or unprovoked venous thromboembolism (Alfonso Iorio. Arch Intern Med 2010; 170:1710-1716). For this reason, patients with venous thromboembolism provoked by a transient risk factor are usually treated with anticoagulant agents for 3 months (Alfonso Iorio. Arch Intern Med 2010; 170:1710-1716), whereas patients with venous thromboembolism that was not associated with a transient risk factor are often treated long-term (Alfonso Iorio. Arch Intern Med 2010; 170:1710-1716). The cumulative risk of recurrence at one, five, and 10 years is 15, 41, and 53 percent, respectively, in patients with an idiopathic venous thromboembolism, compared with 7, 16, and 23 percent in patients with a provoked event (Galioto N J. Am Fam Physician 2011; 83:293-300).

Although it is widely accepted that the risk of recurrence in patients with venous thromboembolism provoked by a transient risk factor is low enough to justify stopping anticoagulant therapy after 3 months, this recurrence risk is not well quantified. Furthermore, the risk of recurrence may not be the same in all patients with venous thromboembolism provoked by a transient risk factor.

D) Risk for Arterial Thrombosis

Arterial thrombosis is a common cause of hospital admission, death, and disability in developed countries (and increasingly in developing nations because of global epidemics of smoking, obesity, and diabetes). It usually follows spontaneous rupture of an atherosclerotic plaque, and may:

-   -   Be clinically silent;     -   Contribute to atherosclerotic progression resulting in coronary         stenosis and stable angina, or lower limb artery stenosis and         claudication;     -   Be present as acute ischemia in the heart (acute coronary         syndromes: unstable angina, myocardial infarction), brain         (transient cerebral ischemic attack or stroke), or limb (acute         limb ischemia).

Traditional risk factors (see table 2) remain the most important markers of arterial disease.

TABLE 2 Risk factors Dyslipemia Current smoker Diabetes Hypertension Abdominal obesity Psychosocial factors

The factor V Leiden and prothrombin G20210A mutations show modest but statistically significant associations with coronary heart disease, stroke, and peripheral arterial events, especially in younger persons (age under 55 years) and in women. The relationship of prothrombin G2010A with arterial thrombosis is very modest (OR 1.32; 95% CI 1.03-1.69) (Kim R J. Am Heart J 2003; 146:948-957). The relationship of factor V Leiden mutation and arterial ischemic events is also modest (OR 1.21; 95% CI 0.99-1.49), patients <55 years old were at greater risk for arterial ischemic event (OR 1.37; 95% CI 0.96-1.97) (Kim R J. Am Heart J 2003; 146:948-957).

There is little evidence that other congenital thrombophilias are associated with increased risk of arterial disease.

Need for New Risk Factors

Despite the above mentioned existence of risk factors and diagnostic tools for early diagnosis arterial thrombosis and venous thromboembolism, which includes deep vein thrombosis and pulmonary embolism, are major causes of morbidity and mortality. Recently, results of the first studies of massive genotyping of the genome in relation to different diseases have been published, including ischemic heart disease and some of its risk factors.

Even among high-risk groups it is not possible to identify individuals who will go on to develop thrombosis and/or venous thromboembolism. Therefore, although several strategies exist for both precise the identification of the risk to develop a thromboembolic event, its prevention or the precise diagnosis of a thromboembolic disease, the goal of preventing the clinical burden of thrombosis and/or thromboembolism has not yet been accomplished (Ruppert A. Current Medical Research & Opinion 2010; 26:2465-2473).

Several attempts have been done to use molecular diagnostics to identify subjects at high risk to develop a thrombotic and/or thromboembolic event. Starting from the finding of a common missense mutation in the factor V gene by Bertina (Bertina R M. Nature 1994; 369:64-67), the description of the prothrombin 20210 3′non-coding sequence variant (Poort S R. Blood 1996; 88:3698-3703), and the thermolabile C677T 5,10-methylene-tetrahydrofolate reductase variant. There is also patent document such as EPO0696325B1 describing the use of mutations in coagulation factors. EPO0696325B1 describes the use of mutations in factor V to identify persons at risk to develop thrombotic event. Or the patent document WO05047533A1 describing a method for detecting the presence or absence of a variant nucleotide in at least two SNP sites associated with thrombosis, said SNP sites selected from the group consisting of factor V Leiden G1691A, Prothrombin (Factor II) G20210A, MTHRF C677T, MTHFR A1298C, factor XIII G4377T, and tissue factor plasma inhibitors (TFPI) C536T.

However, none of the attempts tried until now have proved satisfactory efficacy and the initial enthusiasm for test adoption will need to be tempered by formal evidence of clinical benefit deriving from the test.

Accordingly, there is a need for novel markers, including new genetic markers and specific combinations thereof that would successfully and advantageously predict who is at high risk of developing a thromboembolic disease and/or thromboembolic disease complications such as—but not limited to—deep vein thrombosis, pulmonary embolism, acute coronary syndromes (acute myocardial infarction, unstable angina), stroke, transient ischemic attack or stroke in a way that preventive measures could be implemented to keep that risk at the lowest possible level.

There is also a need for novel markers, including new genetic markers and specific combinations thereof that would successfully and advantageously assist the diagnosis of a thromboembolic disease and/or thromboembolic disease complications such as—but not limited to—deep vein thrombosis, pulmonary embolism, acute coronary syndromes (acute myocardial infarction, unstable angina), stroke, transient ischemic attack or stroke in a way that preventive measures could be implemented to keep that risk at the lowest possible level.

SUMMARY OF THE INVENTION

In a first aspect, the invention provides a method which is suitable to solve the limitations of the methods used nowadays to estimate the thromboembolism risk and/or to diagnose the thromboembolic events for a particular subject.

The method provided according to the present invention solves the limitations comprising the steps of determining in a sample isolated from said subject the presence at least of one of following genetic variants Serpin A10 (protein Z inhibitor) Arg67Stop (rs2232698), Serpin C1 (antithrombin) Ala384Ser (Cambridge II), factor XII C46T (rs1801020), factor XIII Val34Leu (rs5985), Factor II (prothrombin) G20210A (rs1799963), factor V Leiden Arg506Gln (rs6025), factor V Cambridge Arg306Thr, factor V Hong Kong Arg306Gly, ABO blood group rs8176719, ABO blood group rs7853989, ABO blood group rs8176743, and ABO blood group rs8176750 is indicative of the risk of suffering a thromboembolic event (fatal or non-fatal acute myocardial infarction, or stroke, or transient ischemic attack, or peripheral arteriopathy or deep vein thrombosis or pulmonary embolism) which is better than the risk assessment done by the methods nowadays in use.

In another aspect, the invention relates to methods for the establishing the probability of an individual of presenting thromboembolic event based on the presence of one or more of the polymorphisms mentioned above in combination with one or more conventional risk factors, wherein the risk is given.

In another aspect, the invention relates to methods for the establishing the probability of an individual of presenting a recurrent thromboembolic event based on the presence of one or more of the polymorphisms mentioned above in combination with one or more conventional risk factors, sociodemographic and clinical characteristics, wherein the risk is given.

In another aspect, the invention relates to methods for the assistance to the diagnosis of a thromboembolic event based on the presence of one or more of the genetic variants mentioned above in combination with one or more conventional risk factors, sociodemographic and clinical characteristics, wherein the probability for diagnosis is given.

In another aspect, the invention relates to methods for the establishing the need for preventive measurements to prevent the development of a thromboembolic event based on the presence of one or more of the polymorphisms mentioned above in combination with one or more conventional risk factors, sociodemographic and clinical characteristics wherein the risk is given.

“Thromboembolic event” in the context of this application should be understood as the alteration of the hemostasis that leads to the development of a blood clot (thrombo) inside a vascular vessel (artery or vein). The thrombo can even obstruct the vascular vessel completely and/or become detached and obstruct another vascular vessel.

“Thromboembolic event” includes among others the following conditions: arterial thrombosis, fatal- and non-fatal myocardial infarction, stroke, transient ischemic attacks, cerebral venous thrombosis, peripheral arteriopathy, deep vein thrombosis and pulmonary embolism.

“Thromboembolic event” in the context of this application is used interchangeably with “thromboembolism”.

“Thromboembolic event” in the context of this application is used interchangeably with “thrombosis”.

“Thromboembolic event” in the context of this application is used interchangeably with “thromboembolic complication”.

“Thrombophilia” in the context of this application should be understood as the disorders of hemostasis that predispose to thrombosis. Included are heritable deficiencies of the natural anticoagulants antithrombin, protein C, and protein S and common mutations in the genes encoding clotting factors and acquired thrombophilias such as antiphospholipid antibodies.

“Thrombo-prophylactic” in the context of this application should be understood as any preventive measure or medication that reduces the likelihood of the formation of blood clots. Thrombo-prophylactic treatment includes administering to a subject in need of such treatment anticoagulant and/or antithrombotic therapy, including heparins (low molecular weight heparins such as enoxaparin, dalteparin, danaparoid, nadroparin, tinzaparin, or ardeparin), vitamin K antagonists (e.g., Warfarin), direct thrombin inhibitors (e.g., argatroban, lepirudin, desirudin, bivalirudin, Dabigatran etexilate or Dabigatran), and/or factor Xa inhibitors (fondaparinux, rivaroxiban, betrixaban, edoxaban, otamixaban, letaxaban, eribaxaban, or apixaban), and as described elsewhere herein. Other thrombo-prophylactic treatment can include mechanical methods such as passive devices (e.g., knee- or thigh-high graduated compression (elastic) stockings), active devices (e.g., external pneumatic compress or intermittent pneumatic compression) devices, or venous foot pumps.

The terms “disease” and “disorder” shall be interpreted in the context of this application interchangeably.

“Mutation” in the context of this application should be understood as the change of the structure of a gene, resulting in a variant form which may be transmitted to subsequent generations, caused by the alteration of single base units in DNA, or the deletion, insertion, or rearrangement of larger sections of genes or chromosomes.

“Genetic variants” in the context of this application refers to genetic differences both within and among populations. There may be multiple variants of any given gene in the human population (alleles), leading to polymorphism.

The terms “polymorphism” and “single nucleotide polymorphism” (SNP) are used herein interchangeably and relate to a nucleotide sequence variation occurring when a single nucleotide in the genome or another shared sequence differs between members of species or between paired chromosomes in an individual. A SNP can also be designated as a mutation with low allele frequency greater than about 1% in a defined population. Single nucleotide polymorphisms according to the present application may fall within coding sequences of genes, non-coding regions of genes or the intronic regions between genes.

The term “sample”, as used herein, refers to any sample from a biological source and includes, without limitation, cell cultures or extracts thereof, biopsied material obtained from a mammal or extracts thereof, and blood, saliva, urine, feces, semen, tears, or other body fluids or extracts thereof.

“Conventional risk factors” in the context of this application should be understood as those described in tables 1 and 2.

“Sociodemographic and clinical characteristics” in the context of this application should be understood as age, gender, diabetes mellitus, smoking, family history of thromboembolic event, pregnancy, and body mass index.

In a further aspect, the invention relates to a computer program or a computer-readable media containing means for carrying out any of the methods of the invention.

In yet a further aspect, the invention relates to a kit comprising reagents for detecting the genetic variants Serpin A10 (protein Z inhibitor) Arg67Stop (rs2232698), Serpin C1 (antithrombin) Ala384Ser (Cambridge II), factor XII C46T (rs1801020), factor XIII Val34Leu (rs5985), Factor II (prothrombin) G20210A (rs1799963), factor V Leiden Arg506Gln (rs6025), factor V Cambridge Arg306Thr, factor V Hong Kong Arg306Gly, ABO blood group rs8176719, ABO blood group rs7853989, ABO blood group rs8176743, and ABO blood group rs8176750.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1A shows that in the absence of a DNA target, a HairLoop™ is held in the closed state.

FIG. 1B shows that when a target binds perfectly (no mismatch) to its HairLoop™, the greater stability of the probe-target duplex forces the stem to unwind, resulting in an opening of the HairLoop™.

DETAILED DESCRIPTION OF THE INVENTION

The authors of the present invention have solved the problems identified above in the methods in use nowadays for the calculation of the risk in a subject to develop a thromboembolic event, as this term has been defined above.

The authors of the present invention have identified a series of genetic variants which are associated with a risk of presenting a thromboembolic event. These genetic variants show predictive and diagnostic value.

Method for Solving the Limitations of the Methods to the Prediction of the Risk to Develop a Thromboembolic Event or for the Diagnosis of a Thromboembolic Event.

The present application solves the above-described limitation of the methods used nowadays to calculate the thromboembolic event risk and/or to diagnosis a thromboembolic event. A particular combination (as described above) of genetic markers is used, selected and evaluated by the inventors after a complex and genuine analysis of thousands of possible markers. Of the different possibilities to construct a genetic risk score (GRS), the inventors have selected a particular one because it provided the best possible results. To calculate the genetic risk punctuation, the accumulated number of risk allele risk from those SNPs listed in table 3 that are present in each individual is considered. For each of the variants studied, every individual can have 0, 1 or 2 alleles of risk. On having calculated the summatory of risk alleles accumulated in the different set of the selected variants (n=12), for each individual a score that could go from 0 to 24 was given. The inventors have generated new algorithms for thromboembolic risk estimation.

The list of polymorphisms which are used in this method of the present invention is given in Table 3.

TABLE 3 SEQ dbSNP Position Position ID DRFAT dbSNP accession in GRCh37 in GRCh38 NO: GENE Variant PROTEIN NAME sequence Allele number Chrom Chrom Strand Chrom  1 FXII (-4T > C) NP FXII GGACGGA T Risk rs1801020  5 176836532 (-) 177409531 GCCATGA  2 GGACGGA C Non- GCCATGA risk  3 AB0 261delG Val87_ NP CTCGTGGT Risk rs8176719  9 136132908 (-) 133257521 Thr88 = fs G ACCCCT  4 CTCGTGGT Non- _ACCCCT risk  5 AB0  526C > G R176G NP TAGGCGC Risk rs7853989  9 136131592 (+) 133256205 C CACCTC  6 TAGGCGC Non- G CACCTC risk  7 AB0 703G > A G235S NP TGCACCCC Risk rs8176743  9 136131415 (-) 133256028 G GCTTCTA C  8 TGCACCCC Non- A GCTTCTA risk C  9 AB0 1059delC P354R NP TCCGGAA C Risk rs8176750  9 136131059 (-) 133255672 CCGTGAGC 10 TCCGGAA_ Non- CCGTGAGC risk 11 SERPIN  262C > T R88X SERPIN CCTGCTG T Risk rs2232698 14  94756669 (-)  94290332 A10 A10 GAAAGATC T 12 CCTGCTG C Non- GAAAGATC risk T 13 SERPIN 1246G > T A416S Cambridge GAAGCAGC Risk rs121909548  1 173873176 (-) 173904038 C1 II T T CAAGTA C 14 GAAGCAGC Non- T G CAAGTA risk C 15 FV 1601G > A R534Q FV CCTGGACA Risk rs6025  1 169519049 (-) 169549811 Leiden GGC A AGGA ATA 16 CCTGGACA Non- GGC G AGGA risk ATA 17 FV 1001G > C R334T FV  GAAAACCA Risk rs118203906  1 169524537 (-) 169555299 Cambridge C GAATCTT AAG 18 GAAAACCA Non- G GAATCTT risk AAG 19 FV 1000A > G R334G FV Hong AAGAAAAC Risk rs118203905  1 169524538 (-) 169555300 Kong C G GGAATC TTA 20 AAGAAAAC Non- C A GGAATC risk TTA 23 FXIII  103G > T V35L VAL35L GCTTCAGG Risk rs5985  6   6318795 (-)   6318562 GC T TGGTG CC 24 GCTTCAGG Non- GC G TGGTG risk CC 25 FII * 97G > A NP Phro- ACTCTCAG Risk rs1799963 11  46761055 (+)  46739505 thrombin Y A AGCCTC AA 26 ACTCTCAG Non- Y G AGCCTC risk AA

When prediction models are used, as for instance, for making treatment decisions, predictive risks may be categorized by using risk cutoff thresholds.

Those skilled in the art will readily recognize that the analysis of the nucleotides present according to the method of the invention in an individual's nucleic acid can be done by any method or technique capable of determining nucleotides present in a polymorphic site. As it is obvious in the art, the nucleotides present in the polymorphic markers can be determined from either nucleic acid strand or from both strands.

Once a biological sample from a subject has been obtained (e.g., a bodily fluid, such as urine, saliva, plasma, serum, or a tissue sample, such as a buccal tissue sample or a buccal cell) detection of a sequence variation or allelic variant SNP is typically undertaken. Virtually any method known to the skilled artisan can be employed. Perhaps the most direct method is to actually determine the sequence of either genomic DNA or cDNA and compare these sequences to the known alleles SNPs of the gene. This can be a fairly expensive and time-consuming process. Nevertheless, this technology is quite common and is well known.

Any of a variety of methods that exist for detecting sequence variations may be used in the methods of the invention. The particular method used is not important in the estimation of cardiovascular risk or treatment selection.

Other possible commercially available methods exist for the high throughput SNP identification not using direct sequencing technologies, for example, Illumina's Veracode Technology, Taqman® SNP Genotyping Chemistry and KASPar SNP genotyping Chemistry.

A variation on the direct sequence determination method is the Gene Chip™ method available from Affymetrix. Alternatively, robust and less expensive ways of detecting DNA sequence variation are also commercially available. For example, Perkin Elmer adapted its TAQman Assay™ to detect sequence variation. Orchid BioSciences has a method called SNP-IT™ (SNP-Identification Technology) that uses primer extension with labeled nucleotide analogs to determine which nucleotide occurs at the position immediately 3′ of an oligonucleotide probe, the extended base is then identified using direct fluorescence, an indirect colorimetric assay, mass spectrometry, or fluorescence polarization. Sequenom uses a hybridization capture technology plus MALDI-TOF (Matrix Assisted Laser Desorption/Ionization—Time-of-Flight mass spectrometry) to detect SNP genotypes with their MassARRAY™ system. Promega provides the READIT™ SNP/Genotyping System (U.S. Pat. No. 6,159,693). In this method, DNA or RNA probes are hybridized to target nucleic acid sequences. Probes that are complementary to the target sequence at each base are depolymerized with a proprietary mixture of enzymes, while probes which differ from the target at the interrogation position remain intact. The method uses pyrophosphorylation chemistry in combination with luciferase detection to provide a highly sensitive and adaptable SNP scoring system. Third Wave Technologies has the Invader OS™ method that uses proprietary Cleavaseg enzymes, which recognize and cut only the specific structure formed during the Invader process. Invader OS relies on linear amplification of the signal generated by the Invader process, rather than on exponential amplification of the target. The Invader OS assay does not utilize PCR in any part of the assay. In addition, there are a number of forensic DNA testing labs and many research labs that use gene-specific PCR, followed by restriction endonuclease digestion and gel electrophoresis (or other size separation technology) to detect restriction fragment length polymorphisms (RFLPs).

In various embodiments of any of the above aspects, the presence or absence of the SNPs is identified by amplifying or failing to amplify an amplification product from the sample. Polynucleotide amplifications are typically template-dependent. Such amplifications generally rely on the existence of a template strand to make additional copies of the template. Primers are short nucleic acids that are capable of priming the synthesis of a nascent nucleic acid in a template-dependent process, which hybridize to the template strand. Typically, primers are from ten to thirty base pairs in length, but longer sequences can be employed. Primers may be provided in double-stranded and/or single-stranded form, although the single-stranded form generally is preferred. Often, pairs of primers are designed to selectively hybridize to distinct regions of a template nucleic acid, and are contacted with the template DNA under conditions that permit selective hybridization. Depending upon the desired application, high stringency hybridization conditions may be selected that will only allow hybridization to sequences that are completely complementary to the primers. In other embodiments, hybridization may occur under reduced stringency to allow for amplification of nucleic acids containing one or more mismatches with the primer sequences. Once hybridized, the template-primer complex is contacted with one or more enzymes that facilitate template-dependent nucleic acid synthesis. Multiple rounds of amplification, also referred to as “cycles,” are conducted until a sufficient amount of amplification product is produced.

Polymerase Chain Reaction

A number of template dependent processes are available to amplify the oligonucleotide sequences present in a given template sample. One of the best known amplification methods is the polymerase chain reaction. In PCR, pairs of primers that selectively hybridize to nucleic acids are used under conditions that permit selective hybridization. The term “primer”, as used herein, encompasses any nucleic acid that is capable of priming the synthesis of a nascent nucleic acid in a template-dependent process. Primers may be provided in double-stranded or single-stranded form, although the single-stranded form is preferred. Primers are used in any one of a number of template dependent processes to amplify the target gene sequences present in a given template sample. One of the best known amplification methods is PCR, which is described in detail in U.S. Pat. Nos. 4,683,195, 4,683,202 and 4,800,159, each incorporated herein by reference. In PCR, two primer sequences are prepared which are complementary to regions on opposite complementary strands of the target-gene(s) sequence. The primers will hybridize to form a nucleic-acid:primer complex if the target-gene(s) sequence is present in a sample. An excess of deoxyribonucleoside triphosphates is added to a reaction mixture along with a DNA polymerase, e.g. Taq polymerase that facilitates template-dependent nucleic acid synthesis. If the target-gene(s) sequence:primer complex has been formed, the polymerase will cause the primers to be extended along the target-gene(s) sequence by adding on nucleotides. By raising and lowering the temperature of the reaction mixture, the extended primers will dissociate from the target-gene(s) to form reaction products, excess primers will bind to the target-gene(s) and to the reaction products and the process is repeated. These multiple rounds of amplification, referred to as “cycles”, are conducted until a sufficient amount of amplification product is produced.

The amplification product may be digested with a restriction enzyme before analysis. In still other embodiments of any of the above aspects, the presence or absence of the SNP is identified by hybridizing the nucleic acid sample with a primer labeled with a detectable moiety. In other embodiments of any of the above aspects, the detectable moiety is detected in an enzymatic assay, radioassay, immunoassay, or by detecting fluorescence. In other embodiments of any of the above aspects, the primer is labeled with a detectable dye (e.g., SYBR Green I, YO-PRO-I, thiazole orange, Hex, pico green, edans, fluorescein, FAM, or TET). In other embodiments of any of the above aspects, the primers are located on a chip.

In other embodiments of any of the above aspects, the primers for amplification are specific for said SNPs.

Another method for amplification is the ligase chain reaction (“LCR”). LCR differs from PCR because it amplifies the probe molecule rather than producing an amplicon through polymerization of nucleotides. In LCR, two complementary probe pairs are prepared, and in the presence of a target sequence, each pair will bind to opposite complementary strands of the target such that they abut. In the presence of a ligase, the two probe pairs will link to form a single unit. By temperature cycling, as in PCR, bound ligated units dissociate from the target and then serve as “target sequences” for ligation of excess probe pairs. U.S. Pat. No. 4,883,750, incorporated herein by reference, describes a method similar to LCR for binding probe pairs to a target sequence.

Isothermal Amplification

An isothermal amplification method, in which restriction endonucleases and ligases are used to achieve the amplification of target molecules that contain nucleotide 5′-[[alpha]-thio]-triphosphates in one strand of a restriction site also may be useful in the amplification of nucleic acids in the present invention. In one embodiment, loop-mediated isothermal amplification (LAMP) method is used for single nucleotide polymorphism (SNP) typing.

Strand Displacement Amplification

Strand Displacement Amplification (SDA) is another method of carrying out isothermal amplification of nucleic acids which involves multiple rounds of strand displacement and synthesis, i.e., nick translation. A similar method, called Repair Chain Reaction (RCR), involves annealing several probes throughout a region targeted for amplification, followed by a repair reaction in which only two of the four bases are present. The other two bases can be added as biotinylated derivatives for easy detection.

Transcription-Based Amplification

Other nucleic acid amplification procedures include transcription-based amplification systems, including nucleic acid sequence based amplification. In nucleic acid sequence based amplification, the nucleic acids are prepared for amplification by standard phenol/chloroform extraction, heat denaturation of a clinical sample, treatment with lysis buffer and minispin columns for isolation of DNA and RNA or guanidinium chloride extraction of RNA. These amplification techniques involve annealing a primer, which has target specific sequences. Following polymerization, DNA/RNA hybrids are digested with RNase H while double stranded DNA molecules are heat denatured again. In either case the single stranded DNA is made fully double stranded by addition of second target specific primer, followed by polymerization. The double-stranded DNA molecules are then multiply transcribed by a polymerase such as T7 or SP6. In an isothermal cyclic reaction, the RNA's are reverse transcribed into double stranded DNA, and transcribed once against with a polymerase such as T7 or SP6. The resulting products, whether truncated or complete, indicate target specific sequences.

Other amplification methods may be used in accordance with the present invention. In one embodiment, “modified” primers are used in a PCR-like, template and enzyme dependent synthesis. The primers may be modified by labelling with a capture moiety (e.g., biotin) and/or a detector moiety (e.g., enzyme). In the presence of a target sequence, the probe binds and is cleaved catalytically. After cleavage, the target sequence is released intact to be bound by excess probe. Cleavage of the labelled probe signals the presence of the target sequence. In another approach, a nucleic acid amplification process involves cyclically synthesizing single-stranded RNA (“ssRNA”), ssDNA, and double-stranded DNA (dsDNA), which may be used in accordance with the present invention. The ssRNA is a first template for a first primer oligonucleotide, which is elongated by reverse transcriptase (RNA-dependent DNA polymerase). The RNA is then removed from the resulting DNA:RNA duplex by the action of ribonuclease H (RNase H, an RNase specific for RNA in duplex with either DNA or RNA). The resultant ssDNA is a second template for a second primer, which also includes the sequences of an RNA polymerase promoter (exemplified by T7 RNA polymerase) 5′ to its homology to the template. This primer is then extended by DNA polymerase (exemplified by the large “Klenow” fragment of E. coli DNA polymerase I), resulting in a double-stranded DNA (“dsDNA”) molecule, having a sequence identical to that of the original RNA between the primers and having additionally, at one end, a promoter sequence. This promoter sequence can be used by the appropriate RNA polymerase to make many RNA copies of the DNA. These copies can then re-enter the cycle leading to very swift amplification. With proper choice of enzymes, this amplification can be done isothermally without addition of enzymes at each cycle. Because of the cyclical nature of this process, the starting sequence can be chosen to be in the form of either DNA or RNA.

Methods for Nucleic Acid Separation

It may be desirable to separate nucleic acid products from other materials, such as template and excess primer. In one embodiment, amplification products are separated by agarose, agarose-acrylamide or polyacrylamide gel electrophoresis using standard methods (Sambrook et al., 1989, see infra). Separated amplification products may be cut out and eluted from the gel for further manipulation. Using low melting point agarose gels, the separated band may be removed by heating the gel, followed by extraction of the nucleic acid. Separation of nucleic acids may also be effected by chromatographic techniques known in the art. There are many kinds of chromatography which may be used in the practice of the present invention, including adsorption, partition, ion-exchange, hydroxyapatite, molecular sieve, reverse-phase, column, paper, thin-layer, and gas chromatography as well as HPLC. In certain embodiments, the amplification products are visualized. A typical visualization method involves staining of a gel with ethidium bromide and visualization of bands under UV light. Alternatively, if the amplification products are integrally labeled with radio- or fluorometrically-labeled nucleotides, the separated amplification products can be exposed to X-ray film or visualized with light exhibiting the appropriate excitatory spectra.

Alternatively, the presence of the polymorphic positions according to the methods of the invention can be determined by hybridisation or lack of hybridisation with a suitable nucleic acid probe specific for a polymorphic nucleic acid but not with the non-mutated nucleic acid. By “hybridize” is meant a pair to form a double-stranded molecule between complementary polynucleotide sequences, or portions thereof, under various conditions of stringency. For example, stringent salt concentration will ordinarily be less than about 750 mM NaCl and 75 mM trisodium citrate, preferably less than about 500 mM NaCl and 50 mM trisodium citrate, and more preferably less than about 250 mM NaCl and 25 mM trisodium citrate. Low stringency hybridization can be obtained in the absence of organic solvent, e.g., formamide, while high stringency hybridization can be obtained in the presence of at least about 35% formamide, and more preferably at least about 50% formamide. Stringent temperature conditions will ordinarily include temperatures of at least about 30° C., more preferably of at least about 37° C., and most preferably of at least about 42° C. Varying additional parameters, such as hybridization time, the concentration of detergent, e.g., sodium dodecyl sulfate (SDS), and the inclusion or exclusion of carrier DNA, are well known to those skilled in the art. Various levels of stringency are accomplished by combining these various conditions as needed. In a preferred embodiment, hybridization will occur at 30° C. in 750 mM NaCl, 75 mM trisodium citrate, and 1% SDS. In a more preferred embodiment, hybridization will occur at 37° C. in 500 mM NaCl, 50 mM trisodium citrate, 1% SDS, 35% formamide, and 100 [mu]g/ml denatured salmon sperm DNA (ssDNA). In a most preferred embodiment, hybridization will occur at 42° C. in 250 mM NaCl, 25 mM trisodium citrate, 1% SDS, 50% formamide, and 200 [mu]g/ml ssDNA. Useful variations on these conditions will be readily apparent to those skilled in the art.

For most applications, washing steps that follow hybridization will also vary in stringency. Wash stringency conditions can be defined by salt concentration and by temperature. As above, wash stringency can be increased by decreasing salt concentration or by increasing temperature. For example, stringent salt concentration for the wash steps will preferably be less than about 30 mM NaCl and 3 mM trisodium citrate, and most preferably less than about 15 mM NaCl and 1.5 mM trisodium citrate. Stringent temperature conditions for the wash steps will ordinarily include a temperature of at least about 25° C., more preferably of at least about 42° C., and even more preferably of at least about 68° C. In a preferred embodiment, wash steps will occur at 25° C. in 30 mM NaCl, 3 mM trisodium citrate, and 0.1% SDS. In a more preferred embodiment, wash steps will occur at 42° C. in 15 mM NaCl, 1.5 mM trisodium citrate, and 0.1% SDS. In a more preferred embodiment, wash steps will occur at 68° C. in 15 mM NaCl, 1.5 mM trisodium citrate, and 0.1% SDS. Additional variations on these conditions will be readily apparent to those skilled in the art. Hybridization techniques are well known to those skilled in the art and are described, for example, in Benton and Davis (Science 196: 180, 1977); Grunstein and Hogness (Proc. Natl. Acad. Sci., USA 72:3961, 1975); Ausubel et al. (Current Protocols in Molecular Biology, Wiley Interscience, New York, 2001); Berger and Kimmel (Guide to Molecular Cloning Techniques, 1987, Academic Press, New York); and Sambrook et al., Molecular Cloning: A Laboratory Manual, Cold Spring Harbor Laboratory Press, New York, 1989.

Nucleic acid molecules useful for hybridisation in the methods of the invention include any nucleic acid molecule which exhibits substantial identity so as to be able to specifically hybridise with the target nucleic acids. Polynucleotides having “substantial identity” to an endogenous sequence are typically capable of hybridizing with at least one strand of a double-stranded nucleic acid molecule. By “substantially identical” is meant a polypeptide or nucleic acid molecule exhibiting at least 50% identity to a reference amino acid sequence or nucleic acid sequence. Preferably, such a sequence is at least 60%, more preferably 80% or 85%, and more preferably 90%, 95% or even 99% identical at the amino acid level or nucleic acid to the sequence used for comparison. Sequence identity is typically measured using sequence analysis software (for example, Sequence Analysis Software Package of the Genetics Computer Group, University of Wisconsin Biotechnology Center, 1710 University Avenue, Madison, Wis. 53705, BLAST, BESTFIT, GAP, or PILEUP/PRETTYBOX programs). Such software matches identical or similar sequences by assigning degrees of homology to various substitutions, deletions, and/or other modifications. Conservative substitutions typically include substitutions within the following groups: glycine, alanine; valine, isoleucine, leucine; aspartic acid, glutamic acid, asparagine, glutamine; serine, threonine; lysine, arginine; and phenylalanine, tyrosine. In an exemplary approach to determining the degree of identity, a BLAST program may be used, with a probability score between e<“3> and e<“100> indicating a closely related sequence.

A detection system may be used to measure the absence, presence, and amount of hybridization for all of the distinct sequences simultaneously. Preferably, a scanner is used to determine the levels and patterns of fluorescence.

Another method for detecting sequence variations is based on the amplification by PCR of specific human targets and the subsequent detection of their genotype by hybridization to specific Hairloop™ probes spotted on a microarray.

HairLoop™ is a stem-loop, single-stranded DNA molecule consisting of a probe sequence embedded between complementary sequences that form a hairpin stem. The stem is attached to the microarray surface by only one of its strands. In the absence of a DNA target, the HairLoop™ is held in the closed state (FIG. 1A). When the target binds perfectly (no mismatch) to its HairLoop™, the greater stability of the probe-target duplex forces the stem to unwind, resulting in an opening of the HairLoop™ (FIG. 1B). Due to these unique structural and thermodynamic properties, HairLoop™ offer several advantages over linear probes, one of which is their increased specificity differentiating between two DNA target sequences that differ by as little as a single nucleotide.

HairLoop™ act like switches that are normally closed, or “off”. Binding to fluorescent DNA target induces conformational changes that open the structure and as a result after washing, the fluorescence is visible, or “on”.

One HairLoop™ is designed to be specific to one given allele. Thus, assessment of a point mutation for a bi-allelic marker requires two HairLoop™; one for the wild-type allele, and one for the mutant allele. The specific sequences for the detection of the polymorphisms described in table 3 using the HairLoop technology are given in table 4.

TABLE 4 SEQ Sequence Position ID comprising the dbSNP in NO: Variant name polymorphism Allele accession Chrom GRCh38 Strand  1 FXII, (-4T > C) GGACGGA T GCCATGA Risk rs1801020  5 176836532 (-)  2 FXII, 46 C > T GGACGGA C GCCATGA Non-risk (+)  3 AB0, 261delG CTCGTGGT G ACCCCT Risk rs8176719  9 136132908 (-)  4 CTCGTGGT_ACCCCT Non-risk  5 AB0, 526C > G GGAGGTG C GCGCCT Risk rs7853989  9 136131592 (+)  6 GGAGGTG G GCGCCT Non-risk 25 AB0, 526C > G TAGGCGC C CACCTC Risk rs7853989  9 136131592 (+) 26 TAGGCGC G CACCTC Non-risk  7 AB0, 703G > A TGCACCCC G GCTTCTAC Risk rs8176743  9   1361314 (-)  8 TGCACCCC A GCTTCTAC Non-risk (+)  9 AB0, 1059delC TCCGGAA C CCGTGAGC Risk rs8176750  9 136131059 (-) 10 TCCGGAA_CCGTGAGC Non-risk (+) 11 SERPIN A10, 262C > T CCTGCTG T GAAAGATCT Risk rs2232698 14  94756669 (-) 12 SERPIN A10, 728C > T CCTGCTG C GAAAGATCT Non-risk (+) 13 SERPIN C1, 1246G > T CGGTACTTG A AGCTGCTT Risk rs121909548  1 173873176 (-) 14 CGGTACTTG C AGCTGCTT Non-risk 27 SERPIN C1, 1246G > T GAAGCAGCT T CAAGTAC Risk rs121909548  1 173873176 (-) 28 GAAGCAGCT G CAAGTAC Non-risk 15 FV, 1601G > A ATTCCT T GCCTGTCC Risk rs6025  1 169519049 (-) 16 ATTCCT C GCCTGTCC Non-risk 29 FV, 1601G > A CCTGGACAGGC A AGGAATA Risk rs6025  1 169519049 (-) 30 CCTGGACAGGC G AGGAATA Non-risk 17 FV, 1001G > C GAAAACCA C GAATCTTAAG Risk rs118203906  1 169524537 (-) 18 GAAAACCA G GAATCTTAAG Non-risk (+) 19 FV, 1000A > G AAGAAAACC G GGAATCTTA Risk rs118203905  1 169524538 (-) 20 AAGAAAACC A GGAATCTTA Non-risk (+) 21 FXIII, 103G > T GGGCACCA C GCCCTGA Risk rs5985  6   6318795 (-) 22 GGGCACCA A GCCCTGA Non-risk 31 FXIII, 103G > T GCTTCAGGGC T TGGTGCC Risk rs5985  6   6318795 (-) 32 GCTTCAGGGC G TGGTGCC Non-risk 23 FII, * 97G > A TCTCAGC A AGCCTCAAT Risk rs1799963 11  46761055 (+) 24 TCTCAGC G AGCCTCAAT Non-risk 33 FII, * 97G > A ACTCTCAGY A AGCCTCAA Risk rs1799963 11  46761055 (+) 34 ACTCTCAGY G AGCCTCAA Non-risk Method to Establish in a More Appropriate Way the Risk Status.

Another object of the present invention is the development of an algorithm to estimate the risk to develop and/or to being suffering a thromboembolic event. The algorithm is shown as functions 1.

Function 1

Estimating the Risk of Thrombosis.

The individual estimation of the risk of thrombosis is based on a logistic regression model. The aim of this model is to calculate the probability that a person has of presenting venous thrombosis according to his/her genetic, sociodemographic and clinical characteristics. To calculate this probability we use the following equation: Probability(Y=1|x ₁ . . . x _(n))=1/1+exp(β₀+β₁ x ₁+ . . . +β_(n) x _(n+)β_(f·g) x _(f) ·x _(g)+ . . . +β_(h·i) x _(h) ·x _(i)), wherein:

-   -   Probability (Y=1|x₁, . . . , x_(n))=probability of presenting a         thrombosis in a particular individual with concrete and         measurable characteristics in a number of variables 1, . . .         , n. This probability could range between 0 and 1;     -   Exp=exponential natural base;     -   β₀=coefficient that defines the risk (the probability) of         thrombosis non related with the variables 1 to n. This         coefficient can take a value from −∞ to +∞ and is calculated as         the natural logarithm of the incidence of venous thrombosis in         the population;     -   β₁=regression coefficient that expresses the risk (higher or         lower) to present thrombosis associated with the value/presence         of the predictor variable x₁. This coefficient can take a value         from −∞ to +∞;     -   x₁=value taken by the predictor variable x₁ in an individual.         The range of possible values depends on the variable;     -   β_(n)=regression coefficient that expresses the risk (higher or         lower) to present thrombosis associated with the value/presence         of the predictor variable x_(n). This coefficient can take a         value from −∞ to +∞;     -   x_(n)=value taken by the predictor variable x_(n) in an         individual. The range of possible values depends on the         variable.

In addition, the model includes the effect of the combination of some variables in terms of interaction or modification of the effect. That is, the effect size (regression coefficient) of a single variable (x_(f)) can be β_(f) but if this variable is present in combination with another variable (x_(g)) the effect size may vary (increase or decrease) and therefore to consider the effect size of the variable x_(f) we will have to consider not only the β_(f) but also a second regression coefficient β_(f·g) by adding the β_(f) and the β_(f·g). Thus:

-   -   β_(f·g)=regression coefficient that expresses the risk (higher         or lower) to present thrombosis associated with the combined         presence of the predictor variables x_(f) and x_(g). This         coefficient can take a value from −∞ to +∞;     -   x_(f)=value taken by the predictor variable x_(f) in an         individual. The range of possible values depends on the         variable;     -   x_(g)=value taken by the predictor variable x_(f) in an         individual. The range of possible values depends on the         variable;     -   β_(h·i)=regression coefficient that expresses the risk (higher         or lower) to present thrombosis associated with the combined         presence of the predictor variables x_(h) and x_(i). This         coefficient can take a value from −∞ to +∞;     -   x_(h)=value taken by the predictor variable x_(h) in an         individual. The range of possible values depends on the         variable;     -   x_(i)=value taken by the predictor variable x_(i) in an         individual. The range of possible values depends on the         variable;

If the patient does not present any mutation or genetic variant of risk but he/she presents a positive family history of venous thrombosis we will include this variable in the model. The regression coefficient of this variable is 1,185 with a range of possible values from 0.200 to 2.500.

The variables included in the model and the regression coefficients of each of these variables are shown in table 5.

TABLE 5 Regression Regression Risk Regression coefficient coefficient Variable exposure coefficient lower limit upper limit Clinical Age <55 y No 0 55-64 Yes 0.811 0.100 3.000 65-74 Yes 1.409 0.100 3.000 75-84 Yes 1.681 0.100 3.000 >84 Yes 2.534 0.500 6.000 Male Yes 0.336 0.050 1.500 Diabetes Yes 0.351 0.050 1.500 Smoking Yes 0.166 0.050 1.500 Body mass index <25 kg/m² No 0 0 0 25-29.9 kg/m² Yes 0.412 0.050 1.500 ≥30 kg/m² Yes 0.820 0.100 3.000 Use of oral contraceptives Yes 1.131 0.100 3.000 Pregnancy Yes 1.435 0.100 3.000 Family history of thrombosis * Yes 1.185 0.100 3.000 Genetic Factor V Leiden Heterozygote AG 0.993 0.100 3.000 Factor V Leiden Homozygote AA 2.890 0.500 6.000 Prothrombin AG 0.293 0.050 1.500 Serpin 10 TG 1.358 0.100 3.000 Factor XII TC/TT 1.633 0.100 3.000 Factor XIII GT/GG 0.198 0.050 1.500 SerpinC TG 2.277 0.500 6.000 ABO (A1 allele) (see table 6) 0.956 0.100 3.000 Interactions Factor V Leiden · Prothrombin AG · AG 1.114 0.100 3.000 Factor V Leiden · ABO AG · (see table 6) 0.599 0.100 3.000 Factor V Leiden · Oral AG · Yes 0.028 0.005 1.000 Contraceptives Factor V Leiden · Pregnancy AG · Yes 1.191 0.100 3.000 Prothrombin · Oral AG · Yes 0.542 0.100 3.000 Contraceptives Prothrombin · Pregnancy AG · Yes 1.673 0.100 3.000 Prothrombin · BMI ≥30 kg/m² AG · Yes 0.772 0.100 3.000 Oral contraceptives · BMI ≥30 kg/m² Yes · Yes 1.218 0.100 3.000 * Only included in the model if the patient does not present any mutation or genetic variant of risk but he/she presents a positive family history of venous thrombosis.

TABLE 6 Definition of A1 allele The subject carries at least one A1 allele at ABO locus is any of the following combinations is present: Genotypes Combination rs8176719 rs7853989 rs8176743 rs8176750 1 GG CC GG CC 2 GG CC GG CdelC 3 GG CG GA CG or CC 4 GdelG CC GG CG or CC 5 GG CG GG CC

Surprisingly, the combination of SNP markers included in the present invention and set forth in table 3 and using the function described in function 1 have proved to be capable to establishing the risk to develop a thromboembolic disease or event with a higher accuracy than that obtained using the methods nowadays in use or published functions including genetic information.

Surprisingly, the combination of SNP markers included in the present invention and set forth in table 3 and using the function described in function 1 have proved to be capable to assist in the diagnosis of a thromboembolic disease or event with a higher accuracy than that obtained using the methods nowadays in use or published functions including genetic information.

By the use of the functions described, a personalized risk is obtained for the development of thromboembolic event, in particular fatal- and non-fatal-myocardium infarction, stroke, transient ischemic attack, peripheral arteriopathy, deep vein thrombosis, pulmonary embolism or a combination thereof.

Example 1

Introduction. Thromboembolic disease has an important genetic component. In addition to the classic FV Leiden (FVL) and prothrombin G20210A (PT), new genetic variants associated with this pathology have been identified. The aim of this study was to determine whether a set of genetic variants selected by us (genetic profile) improves the ability of FVL and PT to predict the presence of thrombosis.

Methods. We included two studies (thrombosis) and controls: MARTHA (1,150 cases/801 controls) designed to evaluate the association of FVL and PT with other risk factors, and a study in Spanish population: PE (249 casos/248 controls). The genetic profile analyzed was: FVL, PT, ABO, C46T (F12), A384S (SERPINC1), R67X (SERPINA10). The association between genetic variants and thrombosis was calculated using the OR adjusted for age and sex. The predictive ability was calculated using the c statistic (AUC-ROC) and reclassification (NRI, IDI) observed when using the FVL, PT or when using the genetic profile.

Results.

TABLE 7 Association between variants and thrombosis [OR (95%)] and the proportion of FVL and PT carriers compared with carriers of the genetic profile (only cases). Association between variants and thrombosis [OR (95% IC)] and the sensitivity of FVL and PT compared with the sensitivity of all the variants FVL + PT All variants FVL PT A1 C46T A384S R67X MARTHA 2.3 0.9 1.8 0.9 0.9 2.3 (1.8-2.8) (0.7-1.1) (1.2-2.7) (0.6-1.4) (0.2-3.7) (1.2-4.6) Sensitivity 50.4 87.5 PE 7.2 2.8 2.62 3.1 4.1 2.5 (2.8-18.9) (1.2-6.8) (1.8-3.8) (1.1-18.8) (0.5-36.9) (0.8-8.1) Sensitivity 19.7 71.5

TABLE 8 -statistic and reclassification, NRI (net reclassification improvement) and IDI (integrated discrimination improvement) comparing the use of the genetic profile to FVL and PT. C-statistic (95% IC) and reclassification, NRI [Net Reclassification Improvement] (95% IC) and IDI [Integrated Discrimination Improvement] (95% IC) of FVL + PT compared with all the variants C-Statistic NRI IDI MARTHA PE MARTHA PE MARTHA PE FVL + PT 0.54 0.58 Ref. Ref. Ref. Ref. (0.51-0.57) (0.56-0.62) All variants 0.58 0.69 5.3 23.4 1 5.9 (0.55-0.61) (0.64-0.73) (−1.1-11.6) (11.1-35.6) (0.3-1.8) (3.71-7.88) p-value <0.001 <0.001 >0.05 <0.001 <0.5 <0.001

Discussion. We demonstrate that the selected genetic profile significantly improves the prediction of the risk of thrombosis, identifying a genetic risk of presenting a thromboembolic event in 51.6% of people who had a thromboembolic event and through analysis of FVL and PT were not at risk genetic. The genetic profile in clinical practice will improve the diagnosis, prevention and treatment of thromboembolic disease. 

The invention claimed is:
 1. A method comprising: a) determining in a sample isolated from a human subject the presence or absence of each of the following polymorphisms: Serpin A10 (protein Z inhibitor) Arg67Stop (rs2232698), Serpin CI (antithrombin) Ala384Ser (Cambridge II), factor XIIC46T (rsl801020), factor XIII Va134Leu (rs5985), Factor II (prothrombin) G20210A (rsl799963), factor V Leiden Arg506Gln (rs6025), factor V Cambridge Arg306Thr, factor V Hong Kong Arg306Gly, ABO blood group rs8176719, ABO blood group rs7853989, ABO blood group rs8176743, and ABO blood group rs8176750; b) detecting the presence of at least one of said polymorphisms in said subject; c) identifying the subject with at least one of said polymorphisms for treatment and d) administering to the subject identified in step c) an anticoagulant and/or antithrombotic therapy.
 2. The method according to claim 1, further comprising determining one or more cardiovascular disease or disorder risk factors in the human subject selected from the group consisting of age, race, sex, body mass index, smoking status, systolic blood pressure, diastolic blood pressure, hospitalization, plaster cast immobilization, surgery, trauma, oral contraceptives or hormone therapy, pregnancy, prolonged travel (≥2 hours), collagen vascular diseases, heart failure, malignancy, medications, myelo proliferative disorders, nephrotic syndrome, recurrent pregnancy loss, abdominal obesity, diabetes mellitus, low density lipoprotein (LDL)-cholesterol level, high density lipoprotein (HDL)-cholesterol level, cholesterol level, triglyceride levels, family history of thromboembolic event, pregnancy, and body mass index.
 3. The method according to claim 1, wherein the sample is an oral tissue sample, scraping, or wash or a biological fluid sample.
 4. The method of claim 3, wherein the sample is saliva, urine or blood.
 5. The method according to claim 1, wherein the presence or absence of the polymorphism is identified by amplifying or failing to amplify an amplification product from the sample.
 6. The method according to claim 5, wherein the amplification product is digested with a restriction enzyme before analysis.
 7. The method according to claim 1, wherein the polymorphism is identified by hybridizing nucleic acid in the sample with a primer label which is a detectable moiety.
 8. The method according to claim 1, where in the subject has a thromboembolic disease selected from the group of fatal or non-fatal myocardial infarction, stroke, transient ischemic attacks, peripheral arteriopathy, vein thrombosis, deep vein thrombosis, pulmonary embolism or a combination thereof.
 9. The method according to claim 1, wherein the anticoagulant and/or antithrombotic therapy comprises administering to the subject heparins, a vitamin K antagonist, direct thrombin inhibitor, factor Xa inhibitor.
 10. The method according to claim 9, wherein the vitamin K antagonist is Warfarin.
 11. The method according to claim 9, wherein the direct thrombin inhibitor is argatroban, lepirudin, desirudin, bivalirudin, Dabigatran etexilate or Dabigatran.
 12. The method according to claim 9, wherein the factor Xa inhibitor is fondaparinux, rivaroxiban, betrixaban, edoxaban, otamixaban, letaxaban, eribaxaban, or apixaban.
 13. The method of claim 1, wherein the presence of at least two polymorphisms are detected in the subject. 